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I. 


Introduction 


The  Naval  Postgraduate  School  has  started  a  project  to 
design  and  build  a  mini-satellite  (ORION)  .  This  mini-satellite 
bus  will  support  a  small  payload  for  various  space  science, 
technology  and  military  missions.  The  first  specific  payload  is 
not  determined,  yet,  and  therefore,  its  major  design  parameters, 
such  as,  launch  vehicle,  attitude  stabilization  scheme,  are  not 
known.  A  feasibility  study  and  a  preliminary  design  oJ  the  mini¬ 
satellite  were  performed  by  Boyd  (Ref.  1)  .  In  this  preliminary 
design  the  launch  of  the  satellite  from  the  GAS  (Get-Away- 
Specials)  canister  of  a  Space  Shuttle  and  single  spin 
stabilization  were  assumed 

This  preliminary  design,  as  defined  in  Ref.  1,  is 
investigated  from  the  structural  dynamics  and  attitude  dynamics 
points  of  view.  These  include  the  survivability  of  the  satellite 
under  severe  shock  and  vibration  environment  during  launch 
periods,  and  the  stability  of  the  single  spin  attitude 
stabilization  of  the  satellite.  A  rigorous  investigation  of  the 
survivability  satellite  under  shock  and  vibration  environment  is 
not  possible  in  this  preliminary  design  stage  since  the 
acceptible  acceleration  level  of  each  electronic  and  mechanical 
components  and  detail  design  of  connections  and  fasteners  are  not 
defined.  However,  it  seems  that  the  present  design  of  bolts  and 
nuts  connection  and  no  damping  treatment  (Figure  1)  contains  high 
failure  probability  areas  considering  the  wide  input  frequency 
range  (20-2000  Hz)  of  random  vibration  in  GAS  Canisters  (Figure 
2)  . 

The  main  body  of  the  NPS  mini-satellite  is  long  cylindrical 
shape  (See  Figure  1) .  The  principal  moment  of  inertia  about  the 
axial  axis  of  symmetry  is  smaller  than  the  principal  moments  of 
inertia  about  the  radial  axes  of  the  main  body  without  the  four 
radial  booms.  It  is  well  known  that  the  attitude  of  a  rigid  body 
spinning  about  the  axis  of  the  minimum  principal  moment  of 
inertia  is  unstable.  Therefore,  in  the  preliminary  design  of  the 
NPS  satellite,  the  lengths  of  booms  and  the  masse  of  the 
magnetometer  at  the  end  of  each  boom  are  such  that  the  total 
moment  of  inertia  about  the  axial  axis  is  the  maximum,  assuming 
the  four  radial  booms  are  rigid.  However,  long  booms  are 
flexible  and  the  flexibility  of  booms  affects  the  attitude 
stability  of  a  spinning  body  negatively.  Therefore,  each 
specific  design  of  a  spin  stabilized  satellite  with  flexible 
parts  should  be  investigated  to  determine  the  stability 
conditions  between  design  parameters. 

Dynamic  equations  governing  3-dimensional  motions  of  a  body 
(especially  a  flexible  body)  in  space  are  highly  nonlinear  and 
complicated.  Therefore,  a  free  or  forced  motion  and  the 
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stability  of  the  motion  of  a  body  in  3-dimensional  space  can  be 
best  ipvestigated  by  the  computer  simulation  for  each  specific 
case  using  a  highly  capable  computer  code  such  as  DISCOS  (Ref. 
2),  or  TREETOPS  (Ref.  3).  In  this  preliminary  design  stage,  and 
without  an  available  code,  the  stability  of  the  simple  spinning 
motion  of  the  NPS  mini-satellite  about  the  axial  axis  of  symmetry 
is  investigated  analytically  in  Chapter  II.  The  analytical  study 
of  the  spin  stability  of  free  motions  of  elastic  disspative 
systems  was  performed  previously  by  Teixeira-Filho  and  Kane  (Ref. 
4).  The  results  of  this  study  were  applied  to  the  stability 
studies  of  a  spinning  rigid  satellite  that  carries  four 
elastically  mounted  rigid  antennas  (Ref.  5)  and  of  a  spacecraft 
consisting  of  a  rigid  body  and  four  elastic  booms  using  finite 
number  of  beam  elements  (Ref.  6)  .  In  the  present  investigation 
of  the  stability  of  free  spinning  motion  of  the  NPS  mini¬ 
satellite  the  basic  results  of  the  Ref.  4  are  applied,  using 
natural  vibration  modes  of  the  booms  as  the  admissble  functions 
of  internal  coordinates. 


II.  Analysis 

II-l  system  Description 

The  NPS  mini-satellite  S  to  be  analyzed  consists  of  a  rigid 
body  B  of  cylindrical  shape  and  four  flexible  booms  C,  D,  E,  F, 
as  shown  in  Figure  3.  B*  is  the  mass  center  of  B  and  X^, 

are  the  principal  axes  of  inertia  of  B.  When  undeformed,  boom  C 
and  D  lie  on  the  opposite  side  of  B  in  the  X^-X2  plane  and  are 

parallel  to  X2  while  booms  E  and  F  lie  on  opposite  sides  of  B  in 

the  Xj^-X^  plane  and  are  parallel  to  X^ .  The  four  booms  are 

assumed  to  be  rigidly  attached  on  the  cyclindrical  surface  of 
radius  R  of  body  B  in  the  radial  direction  and  are  on  the  X^-X^ 

plane.  Each  boom  has  length  L  and  rectangle-tubular  cross 
section  of  flexural  rigidities  EI^^  about  axis  X^  and  EI^  about 

axis  X^  (for  C,D)  or  about  X2  (for  E,  F) . 

B  has  mass  M  and  the  mass  distribution  in  B  is  assumed  to  be 
axisymmetric  about  X^.  Therefore  among  the  principal  moments  of 

inertia  B^,  B^ ,  B^  of  B  for  the  mass  center  B*  associated  with 

X^ ,  respectively,  B2  and  B^  are  the  same.  Each  boom  has 


X^,  X2, 


mass  and  the  mass  of  a  magnetometer  attached  at  the  tip  of 
each  boom  is  m. 


The  orientation  of  B,  or  the  X^, 


in  a  Newtonian 


reference  frame  N  can  be  specified  in  terms  of  three  external 
coordinates,  0  02,  9  2*  ®3  Euler  angles  such 

that  Xj^,  X2/  X^  coordinates  system  is  brought  into  an  orientation 

in  N  by  successive  rotations  0  3  about  X^^,  X2,  X3  from  the 


N, ,  N 


2' 


N3  direction. 


The  deformed  configuration  of  the  satellite  is  described 
using  4n  internal  coordinates  q^,  q2,  _ q^^^,  in  addition  to  the 

external  coordinates.  The  general  deformation  of  a  boom  C  is 
described  using  the  n  internal  coordinates. 

n 

Wc(x,t)  =  ^<l).(x)q.(t)  (1) 

i=l 

where  W  (x,t)  is  the  deflection  of  a  boom  cross-section  at 
c 

position  X  from  the  boom  root  in  X^^  X2  X3  coordinates,  and 

i  =  l,2,...n  are  admissible  functions  chosen  to  be  the  same  as 
the  nonnal  modes  of  the  boom  as  a  cantilever  beam  with  an  end 
mass  in  this  analysis.  The  general  deformation  of  booms  D,  E,  F, 
are  described  using  the  same  form  of  equation  (1)  with  q^,  i  = 

n+l,...,2n  for  boom  D,  with  q^.,  i  =  2n+i,...,3n  for  boom  E,  with 

g^,  i  =  3n+l,...,4n  for  boom  F,  respectively. 


It  is  assumed  that  no  external  forces  and  moments  are 
applied.  Due  to  the  internal  forces  of  interaction  between  the 
particles  of  systems,  it  is  presumed  that  energy  dissipation  can 
be  described  by  functions  D^{q,q)  which  are  functions  of 

‘^l'"''‘^4n'  ^l'*’*'^4n  such  that 

Dj.(q,0)  =0,  r  =  1,2,  .  .  .  ,4n  (2) 

Dj.(q/q)  kj.  -  0  (3) 


where 

q^ ,  .  .  . 


the  equality  sign  of  equation 
q^^  all  vanish. 


(3) 


holds  only  when 


•  •  • 


1-2  Stability  Conditioi.s 


In  torque-free  motion,  the  satellite  system  S  can  move  such 
that  the  angular  velocity  of  B  in  the  Newtonian  reference  frame  N 
has  a  constant  magnitude  q  and  is  parallel  to  which  remains 

fixed  in  N  while  the  booms  remain  normal  to  This  motion  is 

defined  as  a  simple  spin. 

A  simple  spin  is  called  stable  if  the  system  resumes  this 
motion,  or  acquires  a  motion  closely  resembling  it,  subsequent  to 
every  sufficiently  small  disturbance  that  leaves  the  angular 
momentum  unaltered.  If  the  system  returns  to  the  original  simple 
spin  state  after  a  disturbance  as  time  t  approaches  infinity, 
then  the  simple  spin  is  said  to  be  asymptotically  stable. 

Stability  conditions  can  be  established  by  stating 
conditions  under  which  the  function  F(e ,q,q)  defined  as 

F  =  2(E-E^)  (4) 

is  a  Liapunov  function  (Ref.  7),  where  E  is  the  total  energy  of 
the  system  and  E^  is  the  total  energy  of  the  system  during  the 

simple  spin.  It  was  shown  in  Ref.  4  that  for  a  simple  spin  of  an 
axisymmetric  system. 


1° 

F  =  -5^  {^  (  0^+  9,  )  )  +  Z°,„  q,q,  +  M^q^q,  +  03(6,  q,  q ) 


(5) 


where  H  is  the  magnitude  of  angular  momentum,  I  and  I.  are  the 

s  ^ 

principal  moment  of  inertia  of  the  system  about  the  axial  axis 
and  about  the  transverse  axis,  respectively,  O2(0,q,q)  indicates 

terms  of  third  or  higher  degree  in  its  arguments,  and  superscript 
o  means  its  value  at  the  simple  spin  state. 

is  an  element  of  an  inertia  like  matrix  defined  as 
rs 


I 


iC 

$ 


a 


where  the  symbol  <  >  denote  summation  over  all  particles  of  the 
system  S,  p  is  the  position  vector  of  a  particle  relative  to  the 
mass  center  of  F,  p,^  indicates  its  derivative  about  q^,  m  is  the 

-1  • 

mass  of  the  particle,  and  I  is  the  inverse  of  the  inertia 
dyadic  I. 

Z,  is  a  second  derivative  of  an  energylike  function  Z(q) 
rs  \ 

with  respect  to  q^  and  q^,  which  is  defined  as 

Z(q)  =  V(q)+iH^r‘(q)  (7) 

where  V(q)  is  the  potential  energy  function  of  S. 

From  the  Liapunov  stability  theory,  the  stability 
conditions  are 


F  <  0 


with  equality  sign  holds  only  when  vanish, 


is  positive  definite. 


The  first  condition,  equation  (8)  is  satisfied  from  the 

energy  dissipation,  equations  (2)  and  (3) .  The  second  condition, 

equation  (9)  ,  is  the  well  known  "Maximum  Axis  Rule"  of  the 

stability  condition  of  a  rigid  body,  and  can  be  called  the 

"external  stability  condition"  of  a  deformable  body  system.  This 

condition  can  be  satisfied  when  the  boom  length  L  is  longer  than 

the  critical  length  L  .  ,  which  is  about  60  inches  for  the  NFS 

min 

mini-satellite  (ORION) . 

The  third  stability  condition,  equation  (10),  is  due  to  the 


flexible  deformation  of  booms  and  can  be  called  the  "internal 
stability  condition".  It  is  shown  in  Ref.  4  that  for  a  stress- 
free  spin  about  axis  X^,  that  is, 


or 


o 

II 

O 

> 

;  r  =  1,2, . 

.  4n 

(11) 

r  =  1,  2,.... 

.  4n 

(12) 

from  the  associated  eqations  governing  the  internal  coordinates, 
can  be  expressed  as  follows  with  H  replaced  oy  ^ 


,  =v",  (f  f  +f  f  )  +  lf  };  r=l,2,...4n 

’rs  ’rs  0  0  12,r  12,s  13,r  13, 2  H.J's 


(13) 


r- 1 

S  I 


The  simple  spin  motion  of  the  NFS  satelite  system  S  is  a 
stress-free  spin,  since  it  satisfies  the  condition,  equation  (11) 
or  equation  (12),  by  choosing  the  internal  coordinates  as  in 
equation  (1)  .  Moreover,  by  this  choice  of  internal  coordinates 
using  the  normal  vibration  modes  of  the  cantilevered  boom,  the 

components  of  Z?  in  equation  (13)  can  be  easily  calculated, 

^  S 

using  the  orthonormal  relations  between  normal  modes  of  the  boom: 


m(t)/L)<t)^(L) 


+  J<t)^(x)(>^(x)  1,1  dx  = 
L 


(14) 


jEI(J)^(x)<t)^(x)  dx  = 
L 


1 

I 


(15) 


where  5 


I 

I 

’v 

I 


I 

i 

>$1 

I 

•'ll 

iS 


where  5  is  the  Kronecker  delta  of  r,  s  =  l,2,...,n,  and  m 
and  r  *  l,2,...,n,  are  modal  masses  and  modal  stiffnesses  of 

the  vibration  modes,  y  is  the  mass  per  unit  boom  length. 


V"  =k  5  +K^ 

IS  n  IS  G 


=  {m+(R+L-x))  (l)^(x)  <))’ (X)  , 


where  K_  is  the  geometric  stiffness  of  the  boom  due  to  the 
centrifugal  force  of  the  spin,  and  Wj.  is  the  natural  frequency  of 
the  r  mode. 


IXTl.j 

^2.1  ^  dx  +  L  <|)^(L)  m 


r  =  vibration  modes  in  X^-  plane  for  booms  C  and  D 


IVT-A- 

^3,S  =  n  dx  +  L  (1)^(L)  m 


s  =  vibration  modes  in  X^-  X^  plane  for  booms  E  and  F 


1°  =m  5 

11, IS  rr  IS 


r  =  vibration  modes  in  X^-  X^  plane  for  booms  C,  D,  E, 
and  F 


When  all  ingredients  required  for  the  evaluation  of  are 

available  using  equations  (13)  ,  (16-19)  ,  the  internal  stability 
conditions  can  be  examinded  as  followes: 


I. 


Define  A.  ,  i  =  1,2, 


4n,  as 


(20, a) 


A  ^ 


7° 

^  ’ll 

z° 

^  ’12 

0 

N 

7° 

^  ’22 

z° 

^  ’ll 

z° 

^  ’12 

z“ 

^  ’13 

z° 

^  ’21 

7° 

^  ’22 

z® 

^  ’23 

z° 

^  ’31 

z° 

^  ’32 

z® 

^  ’33 

A  ^ 


and  so  on. 

The  internal  stability  condition  is 


(20, b) 


(20, c) 


Aj  >  0  for  all  i  =  1,2,3, . 4n 


(21) 


II-3  Stability  of  the  NPS  Mini-Satellite  (ORION) 

The  stability  conditions  of  the  simple  spinning  motion  of 
the  NPS  mini-satellite  between  the  angular  velocity  ^  and  the 
length  of  a  boom  L  for  other  parameters  defined  in  the 
preliminary  design  (See  Ref.  1).  The  numerical  value  of  these 
parameters  are  shown  in  Table  l. 


le  1  -  Parameters  of  WPS  mini-satellite 


Radius  of  B 
Mass  of  B 

Axial  moment  of  inertia  of  B 
Transverse  moment  of  inertia  of  B 

Area  of  a  boom  cross-section 

Area  moment  of  inertia  of  a  boom 
about*  X^-  X^  plane 

Area  moment  of  inertia  of  a  boom 
about  X^-  X^  plane 

Young's  modulus  of  the  boom 

Mass  of  a  boom  per  unit  length 

Mass  of  a  magnetometer  at  the  end 
of  a  boom 


R  =  9.5  inch 
M  =  7.542  slug 
Bj  =  2.18  slug-ft 

B2  =  =  6.78  slug-ft 

A  =  0.212  in 

=  0.01668  in 

I2  =  0.06328  in 

E  =  10.3  X  10  Psi 
y  =  0.000415  (Slug/in) 

m  =  0.0621  slug 


For  the  design  value  of  L  =  78.5  inches  the  cantilever  beam 
model  of  a  boom  with  an  end  mass  was  analyzed  using  the 
MSC/NASTRAN  code  for  the  natural  frequencies  and  mode  shapes. 
The  natural  frequencies  are  shown  in  Table  2,  and  the  finite 
element  model  and  the  first  3  mode  shapes  are  shown  in  Figures  4 
to  7. 


Mode 


Natural  Frequency  (Hz) 


Mode  Shape 


1 

2.104 

1-2  plane  bending 

2 

4.097 

2-3  plane  bending 

3 

23.426 

1-2  plane  bending 

4 

45.569 

2-3  plane  bending 

5 

73.656 

1-2  plane  bending 

6 

142.980 

2-3  plane  bending 

When  the  angular  velocity  of  simple  spin  n  is  increased 
from  very  small  value,  it  can  be  seen  from  stability  conditions, 
equations  (20-21)  and  from  equations  (13),  (16-19)  that  the 

stability  condition  will  be  violated  first  when  approaches  the 
first  natural  frequency  since  for  higher  modes  w  ?,  i  = 

2,3,...,  are  much  greater  than  0^ a^d  the  matrix  Z?  become 

5 

diagonally  dominant.  Therefore,  the  stability  conditions  can  be 
examined  simply  by  reducing  the  number  of  internal  coordinates  q^ 

to  the  first  mode  of  each  boom,  q^,  q^,  q^. 


Inserting  the  expressions  of  its  components,  equations  (16) 

-  (19)  into  equation  (13) ,  the  matrix  Z?  becomes  the  following 

iTS 

form: 


= 

rs 


co^  -  aQ^ 
1 


0 

0 


0 

(0^  -  aQ^ 
1 

0 

-PQ2 


0) 


0 


0 


0 

-pQ2 

0 


co^  -  aQ^ 


(22) 
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where 


R+L 


R+L 


a  =  -^{jx<|>j(x)  ^idx  +  L<j)j(L)m}^- {m+^(R+L-x)}{(t)j(x)}^dx  (23) 

IT'I.  R  D 


S  t 


and 


R+L 


P  =  — 1— {fx<t)  (x)ndx  +  L()>  (L)m 
tO  tO  J  ^  1 

1  ■!.  R 


(24) 


s  t 


with  il*j^(x)  is  the  mode  shape  of  the  first  mode. 

The  stability  conditions,  equations  (20)  -  (21) ,  become 


0)^  -  a£l^  >  0 
1 


(25, a) 


(coj  -  >  0 


(25, b) 


(CG?  -  afl2){(o)2  -  ii"*}  >  0 


(25, c) 


Kcaj  -  a£22)2  .  p2Q4j2  >  g 


(25, d) 


Therefore,  the  internal  stability  conditions  is  reduced  to  the 
condition 


Ci>2 


a+P 


(26) 


from  the  condition  (25, c). 


15 


As  the  length  of  a  boon  L  Increases,  the  first  natural 
frequency  decreases  as  follows; 


CO 


The  denominator  in  equations  (23) ,  (24) ,  1°  -  1°  is 

S  w 

.  2  2 

proportional  to  (L  -  L  .  )  and  increases  from  zero  for  L  >  L  .  , 
^  '  min'  min' 

which  is  the  external  stability  condition.  Therefore,  the 

stability  boundary  in  the  stability  diagram  varies  as  shown  in 

Figure  8.  For  the  designed  value  of  L  =  78.5  inches  the 

stability  condition  is 


8.658  (rad/sec)  (28) 


However,  it  was  assumed  that  each  boom  is  a  single  piece  and 
rigidly  connected  to  the  rigid  main  body  B.  Since  each  boom 
consist  of  three  beam  and  the  connection  between  the  boom  and  the 
body  B  is  not  ideally  rigid  the  critical  angular  velocity  ^ 

o 

may  be  much  lower  than  the  value  calculated,  equation  (28) , 
depending  on  the  design  of  joints  of  booms. 


««> 


III.  Conclusions 

The  stability  conditions  of  a  simple  spin  motion  of  the  NFS 
mini-satellite  (ORION)  are  derived  analytically  using  Liapunov 
stability  theory  and  by  choosing  the  natural  vibration  modes  as 
admissible  functions  of  internal  coordinates  of  the  system.  By 
applying  this  conditions  to  the  preliminary  design  of  the 
satellite  it  is  shown  that  the  angular  velocity  of  the  simple 
spin  should  be  limited  to  the  value  of  8.685  rad/sec  for  the  boom 
length  L  =  78.5  inches.  If  the  boom  length  is  changed  the 
critical  angular  velocity  a  also  changes  and  the  stability 

conditions  should  be  re-examined  using  the  results  derived  in 
this  study. 

It  should  be  recognized  that  the  stability  conditions  were 
derived  for  the  simplified  ideal  model  under  the  following 
assumptions;  1)  the  mass  distribution  of  the  cylindrical  main 
body  is  axisymmetric;  2)  the  connections  between  the  booms  and 
the  main  body  and  the  main  body  itself  are  rigid;  3)  each  boom 
acts  as  a  single  piece  and  the  locking  mechanism  of  the  boom 
should  be  ideal. 

Due  to  the  highly  nonlinear  behavior  of  a  three  dimensional 
motion  of  a  flexible  body,  the  stability  analysis  is  valid  for 
only  small  disturbances.  Therefore,  stability  of  dynamic  motion 
of  a  satellite  for  general  environments  should  be  investigated  by 
dynamic  simulation  using  computers  and  proper  dynamic  models. 
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Figure  8.  Stability  Diagram  of  The  NFS  Mini-Satellite. 
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